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ABSTRACT 


-  In  Part  I  of  this  report  the  optimum  MTI  receiver  was  derived  and 
analyzed  for  the  case  in  which  the  radar  pulses  were  emitted  from  the  trans¬ 
mitter  equally  spaced  in  time.  For  typical  long  range  ATC  surveillance  radars, 
aliasing  of  the  target  and  clutter  spectra  results  in  detection  blind  speeds 
at  multiples  of  approximately  70  knots.  It  is  well  known  operationally  that 
these  blind  speeds  can  be  eliminated  by  staggering  the  transmitter  PRF. 
Heretofore,  there  has  been  no  thorough  theoretical  analysis  of  the  effect  of 
staggered  PRF  on  the  spectral  distribution  of  the  target  and  clutter  signals. 

It  is  shown  in  Part  II  that  the  clutter  spectral  density  continues  to  fold  over 
at  the  PRF,  but  that  the  signal  spectrum  becomes  dispersed  in  frequency,  some¬ 
what  like  an  anti-jam  signal.  The  effect  that  this  phenomenon  has  on  the 
performance  of  the  optimum  processor  is  evaluated  in  terms  of  the  signal -to- 
interference  ratio  (SIR)  criterion  that  was  derived  in  Part  I. 

It  is  further  noted  that  even  when  the  target  Ooppler  shifts  are  more 
than  one  PRF  apart,  the  spectra  are  distinguishable,  suggesting  that  unambiguous 
Ooppler  estimation  may  be  possible.  This  concept  is  explored  in  detail  using 
the  MTI  ambiguity  function.  It  is  shown  that  good  SIR  performance  can  be 
obtained  by  choosing  the  stagger  parameters  to  minimize  the  height  of  the 
subsidiary  Doppler  side-lobes.  The  resulting  design  problem  is  noted  to  be 
similar  to  that  of  obtaining  good  antenna  patterns  fo'  arrays  having  non- 
uniformly  spaced  elements. 
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Acting  Chief,  Lincoln  Laboratory  Liaison  Office 

ill 


A  Theory  for  Optimal  Mil  Digital  Signal  Processing 
Part  II:  Signal  Design 

I.  INTRODUCTION  AND  SYNOPSIS 

In  Part  I  of  this  report  [1]  statistical  decision  theoretical  methods 
were  used  to  develop  a  rational  basis  for  comparing  the  performance  of  MTI 
receivers.  The  analysis  has  led  to  the  development  of  a  new  receiver  struc¬ 
ture  that  is  practical  to  implement  using  digital  signal  processing  (DSP) 
techniques  and  achieves  essentially  optimum  performance.  All  of  the  results 
in  Part  I  were  based  on  the  assumption  that  pulses  leave  the  transmitter 
uniformly  spaced  in  time.  For  en-route  L-band  radars  in  which  the  unambigu¬ 
ous  range  must  be  200  n.  mi.,  unambiguous  velocity  measurements  are  not  poss 
ible  because  of  target  spectrum  aliasing  at  the  PRF.  Furthermore,  the  clut¬ 
ter  spectrum  also  folds  over  at  the  PRF  resulting  in  "blind  speeds"  at  which 
the  detection  SNR  of  even  the  optimal  detector  is  degraded  below  practically 
useful  limits.  This  effect  is  demonstrated  in  Figure  1.  In  the  development 
of  classical  MTI  processing  it  has  been  found  from  intuitive  considerations 
that  if  the  transmitter  pulses  are  staggered  in  time,  improved  detection  per 
formance  can  be  obtained  [2],  [3].  However,  there  has  been  no  thorough  theo 
retical  investigation  of  the  exact  effect  that  staggered  PRF's  have  on  the 
underlying  target  and  clutter  models.  The  analysis  developed  in  Part  I  is 
generalized  in  this  report  to  allow  for  the  non-uniformly  spaced  sampling 
pattern.  In  Section  II,  models  are  derived  for  the  sampled-data  target  and 
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clutter  returns  that  result  when  a  ranye  ring  is  sampled  by  repetitive  bursts 
of  M  non-uni formly  spaced  pulses.  The  resulting  model  is  used  in  conjunction 
with  the  decision  theoretical  test  of  Part  I  to  derive  the  optimum  receiver 
structure.  As  in  the  uniformly  sampled  case.the  processor  consists  of  a 
clutter  rejection  filter  and  a  bank  of  matched  filters  that  are  in  a  sense 
matched  to  the  target  signal  over  the  enlarged  unambiguous  velocity  region. 

It  is  shown  that  this  enlarged  complex  of  matched  filters  can  be  realized 
by  making  appropriate  interconnections  of  filters  that  extend  over  only 
the  original  ambiguous  frequency  interval.  Hence,U  may  very  well  be 
practical  to  implement  the  optimum  processor  using  DSP  techniques.  The  Signal 
to-Interference  (SIR)  performance  measure  is  used  to  evaluate  the  performance 
of  the  optimum  detector  and  it  is  shown  that  reasonable  detection  can 
be  achieved  at  velocities  that  previously  could  not  be  seen  by  the  radar. 

In  addition  to  providing  better  detection  performance  over  a  larger  velocity 
interval,  the  optimal  processor  is  capable  of  providing  velocity  estimates 
over  the  larger  velocity  range.  Since  staggering  the  PRF  increases  the 
unambiguous  velocity  interval  at  the  expense  of  a  decrease  in  the  unambiguous 
range  interval,  it  is  clear  that  the  ambiguity  surface  of  the  transmitted 
waveform  is  being  altered.  Therefore, staggering  the  PRF  is  basically 
an  MTI  signal  design  problem  and  hence  is  characterized  by  the  range-velocity 
ambiguity  function.  This  function  is  evaluated  along  the  Doppler  axis 
as  this  represents  the  output  of  the  matched  filters  of  the  optimal  processor 
It  is  shown  that  the  M-  pulse  staggered  waveform  reduces  the  velocity 
ambiguity  at  the  average  PRF. 


II.  INTRODUCTION  TO  MTI  SIGNAL  DESIGN 

The  analysis  presented  in  Part  I  has  led  to  the  development  of  a 
quantitative  technique  for  evaluating  optimal  and  suboptimal  MTI  receivers. 
The  results  show  that  a  considerable  improvement  in  target  detection  capa¬ 
bility  is  possible  using  the  matched  filter  receiver.  The  problem  formula¬ 
tion  and  receiver  synthesis  are  based  on  the  assumption  that  the  sampling 
rate  is  uniform.  In  that  case,  for  the  L-band  ARSR[4],  an  aircraft  moving 
at  600  kts.  induces  a  Doppler  shift  corresponding  to  3000  Hz.  Since  the  PRF 
needed  to  obtain  200  nmi.  unambiguous  range  is  360  pulses/see.,  aliasing  of 
the  target  and  clutter  spectra  will  occur  with  period  360  Hz.  or  72  kts. 
Therefore  if  an  aircraft  is  moving  at  a  velocity  +  n  x  72  kts.  n  =  0,1,2,..., 
the  Signal-To-tnterference-Ratio  (SIR)  will  be  seriously  degraded  due  to  the 
clutter  aliasing.  Furthermore  it  will  be  impossible  to  distinguish  between 
a  target  moving  at  velocity  v  and  another  at  v  +  n  x  72.  Since  staggering 
the  PRF  has  been  found  to  improve  the  detection  capabilities  of  MTI  receivers 
at  the  blind  speeds  [2]  it  is  of  interest  to  determine  the  theoretical  basis 
for  this  improvement  and  to  explore  its  implications  regarding  the  question 
of  velocity  resolution.  Since  the  underlying  statistical  properties  of  the 
data  samples  will  be  affected  by  the  non-uniform  sampling  pattern,  it  is 
necessary  to  re-examine  the  basic  target  and  clutter  models  that  were  derived 
in  Part  I  for  the  uniformly  sampled  system. 
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It  was  shown  in  Section  II  of  Part  I,  that  if  the  aircraft  induced  a 


Doppler  frequency  v  and  was  located  at  azimuth  <j>  =  then  uniformly  spaced 
transmit  pulses  led  to  target  samples  at  a  range  cell  given  by  ^  I -(14),  namely 

j2irvnT 

s(nTp;a)  =y  g(nTp-T)  e  p  (1) 

where  g(t)  is  the  two-way  antenna  voltage  gain  pattern  and  Tp1s  the  uniform 
interpulse  period.  In  the  derivation  of  this  target  model, it  was  assumed 
that  the  transmitted  pulses  were  narrow  compared  to  the  Doppler  period  and 
to  antenna  pattern  variations.  In  other  words, the  physical  sampling  was  done 
by  modulating  a  continuous  phenomenon  by  a  train  of  sampling  pulses.  A  useful 
idealization  is  to  represent  the  sampled  data  sequence  as  the  continuous  time 
function  as  follows: 

CO 

$(t;a)  =  s(t;a)  'S  $(t-nT  )  (2) 

—  y 

n=  -<o 

where 


j2wvt 

${t;a)  BY9(t-t)  e  (3) 

Then  the  Z-Transfor®  of  the  uniformly  sampled  sequence  is  related  to  the 
Fourier  Transform  as  follows: 
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The  notation  I-(14)  refers  to  equation  (14)  in  Part  I. 


•  J  — 


rsr  "^i  -w  «  .-?»-^7jW,  ,^' .  1 1;  u  * 


j2irfTn  * 

Zs(e  p;<*)  -  l  [s(nT  ;a)3 


j2TTfT 


z  =  e 


P  i 


3  F[s(t;a)] 


00 


2$(f "  J- »$*.) 
n=-®  P 


where 


(4) 


S(f ;a)  -  F  [s(t;a)]  (5) 

Equation  (4)  shows  the  foldover  of  the  target  spectrum  every  l/TpHz. 

When  a  two-pulse  staggered  PRF  sampling  pattern  is  used,  samples  of 
the  target  environment  are  taken  at  times  0,  ±  (Tp  -  e),  *  2Tp,  *  (3Tp  -  c), 
*4  Tp,  . as  shown  in  Figure  2.  In  this  case,  the  sequence  of  samples  has 
values 


...  s(o),s(Tp-c).s(2Tp).s(3Tp-c),s(4Tp),  ...  (6) 

These  numbers  correspond  to  sampling  s(t)  at  times  ...  0,  2Tp,  4Tp,  6Tp,  .... 
and  sampling  s{t-e)  at  times  ...  Tp>  3Tp,  5Tp,  7Tp,  •**.  A  continuous  time 
representation  of  the  sampled-data  waveform  for  the  two-pulse  staggered 
algorithm  is  therefore; 
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TRANSMITTER 

OUTPUT 


Fig,  2.  Two-pulse  staggered  sampling  pattern 


00 


oo 


(7) 


s(t;a)  ■  s(t;a)  6(t-2nTp)  +  s(t-e;a)^>  5[t-(2n-l)Tp] 

f|3-CO 


For  the  M-pulse  stagger  sampling  pattern, 

s(t)  Is  sampled  at  0,  MTp,  2MTp,  . . 

s(t-e-j)  is  sampled  at  Tp>  (M+l)  Tp,  (2M+1)  TpS  ... 

s(t-c2)  is  a;..,  'ed  at  2Tp,  (M+2)Tp,  (2M+2)  Tpa  ... 


s(t-eM_1)  is  sampled  at  (M-l)T  ,  (2M-1)  Tp,  (3M-1)  Tp,  ...  (8) 


from  which  it  is  possible  to  deduce  the  following  continuous  time  representa¬ 
tion  of  the  sampled-data  waveform: 

M-l 

s(tia)  •  £  l^t-^ja)  ^  6[t-(nM+m)Tp] j  (9) 

m=0  n-”°o 

2 

The  Fourier  Transform  of  this  function  is  : 


M-l 

S(f;o)  a 

m*0 


-j2irfe 
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J2irfmT 
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00} 


o 

The  asterisk  will  denote  convolution. 
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Use  will  be  made  of  the  following  identity 


r-  j27rfmT 

>,  •  p  t(f  -  rf-j 


«°  M-l 
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i=-°°  k=0 
»  M-l  km 

.  y  y  e'J  ^ 


i=-oo  k=0 
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In  addition,  for  the  target  signal  of  interest 


S(fjfl)  =  Y  Fg(f-v)  e 


~j2rrfT 


where  F^(f )  is  the  Fourier  Transform  of  g(t)  and  y  =  y  e 
(12)  and  (11)  in  (10),  the  target  spectrum  becomes 


J2ttvt 


.  Using 


S(f;o) 


'V  -j2irf(x+em)  f  i  v  ^ 

’  *  2  Fg<f->  *  *  Mr  2  2 


oo  M-l  km 

i  2  e  ^  -  T, 

=-«  k=n  P 


_  L  V  1  V  V  ‘^TT  -  -f^)  (T+em)  k  , 
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Since  the  target  is  sampled  only  at  discrete  instants,  the  delay  parameter  t 
can  be  estimated  only  to  within  an  interpulse  period.  Therefore,  it  can  be 
assumed  that 


t  *  I ( t)  Tp 

where  I { t)  is  s-'me  unknown  integer.  Then  (13)  becomes 

S(fia)  = 


(14) 


'  *j2TTfT  ” 

H  i=-» 


IV  -J2*r  "■j2ir(f+WT-'T")%iF  jj 

le 

V.1Q  e  p  p  Fg(f-V+MT~  ~ 

o 

It 

■w 

.  m=0  p  p  j 

05) 


-j2Trfe_. 

Since  the  term  e  changes  slowly  relative  to  the  width  of  the  function 

Fg(f-v),  then  to  a  good  approximation 


e’^  Fg(f-v>  -  e'"™”  Ft(f  v, 


(16) 


and  (15)  can  be  written  as 


S(f;a) 


*  -JMt" 

V  £ 

v  i=-» 


Ie 

k=0 
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^e'j2TTr  e‘j2irven 


L  m=0 


+  S ■  V 


(17) 
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It  is  appropriate  to  define  the  coefficients 


-j2TT 

ak(x)  *  e 


kl(x) 

nr 


M-l  .0  km  .0 

h ,  ,  i  v  'j2:tr 'j2’’vs 
bk(u)  -m26  e 


m=0 


ck(x,v)  =  ak(x)bk(v) 
k  =  0,  I,  •••,  M-l 


(18a) 

(18b) 

(18c) 


as  this  leads  to  the  following  convenient  expression  for  the  Z-Transforw  of 
the  target: 

j2trfT  ,  -j2irfx 

Zs  (e  p;a)  3  y  e 


T“  V 
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Clutter  Model 

In  Section  II  of  Part  I,  it  was  shown  that  each  clutter  scattering 
center  could  be  treated  as  a  point  target  having  zero  Doppler.  Therefore,  as 
in  I-(15),  the  nth  scatterer  at  azimuth  $n  in  the  particular  ring  of  interest 
generates  the  clutter  signal  return 

c(tk)  .  9(tk  -T„>  (20) 
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where  rn  *  $n/u>s  -  rnAT/2,  yn  -  An  e^0n  and  tn  represents  the  times  the 
samples  are  taken.  As  before,  An  Is  related  to  the  scattering  cross-section 
of  the  scatterer  and  eR  the  carrier  phase  it  introduces.  From  scan-to- 
scan,the  shift  in  transmitter  phase  and  the  jitter  in  the  antenna  rotation 
render  en  and  Af)  random  variables,  but  over  any  one  scan,  (20)  represents 
a  deterministic  signal  return.  Hence,  the  analysis  used  to  derive  the 
Fourier  Transform  of  the  non-uniformly  sampled  target  return  is  directly 
applicable  to  (20).  The  using  (15)  the  transform  is 


C„(f) 


Y„ 

v~  e 
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I 

1»-ga 


M-l 

k=>0 


-j 2tt 

"n — 


M-l 


i  y  -j2irr  'jaT{f^:x)emF  if+k  i  , 

RZe  e  p  p  Fg(fX  v 

rpO  m  - 


i 


(21) 


Equation  (21)  is  derived  from  (15)  rather  than  (19)  because  the  latter 
equation  has  made  use  of  the  approximation  in  (16).  Since  clutter  returns 
can  be  orders  of  magnitude  greater  than  the  signal  returns,  approximations 
cannot  be  made  unless  they  can  be  justified  on  the  basis  of  signal-to-clutter 
ratios.  The  total  clutter  return  is  due  to  a  finite  number  of  scatterers, 
hence 


C<‘l>  •  I  Cn(tk* 
n 


(22) 
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and  the  Fourier  Transform  of  this  aggregate  of  returns  is 


C(f)  -  2,  Cn(f)  (23) 

n 

Therefore,  the  energy  spectral  density  of  the  clutter  measured  over  a  single 
scan  is 


C(f) 
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II  \ 


nl  n2 


(f)  C  *(f) 
n2 


(24) 


This  is  a  random  process  in  the  sense  that  after  each  scan  the  values  of 
and  en  change  in  a  random  fashion.  Then  the  average  power  spectral  density 
of  the  clutter  is 


sc(f) 


C(f) 


1 

*7 


v 

nl 


y  c  (f)  c*  <f) 

is  n%  ru 

m  1  • 


(25) 


where  T,  is  the  scan  time  and  the  bar  denotes  statistical  averaging  over  the 

b 

random  variables  Ar  and  9n> 


Since  the  amplitudes,  phases  and  azimuthal  locations  of  the  scatterers 
are  independent,  each  of  the  random  variables  in  (25)  can  be  averaged 
separately.  Furthertaore,  it  follows  that 


(26) 
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and  since  the  frequency  extent  of  F  (f)  is  narrow  relative  to  a  separation 

w 

k/MTp, 


k-i  i  ■]  it  k«  i  o  » 

Fg<f  *  Hr  ■  r>  Fg  (f  +  r  -  r> » 

y  p  p  y  p  p 


,  k,  1, 

Fs!f  +  fC  -  r> 

P  P 


\  »|y  ^4 

ki  k2  M  h 


(27) 


Substituting  (21)  in  (25)  and  using  (26)  and  (27),  it  follows  that 


Sc(f» 


p  i=- 


M-l  M-l  9  km 
k=0  m=0 


P  P 


(28) 


where 


A  (f)  a  e 
m  ' 


“j2irfem  p 


(f) 


(29) 


In  (29),  l/em  is  generally  much  greater  than  the  frequency  extent  of  the 
clutter  and  It  is  reasonable  to  assume  that 

-j2T.fr 

e  Fg(f)  “  Fg(f'  (30) 

Since  the  clutter  signals  can  be  many  orders  of  magnitude  greater  than  the 
signal,  this  approximation  must  be  undertaken  with  care  in  each  application. 
An  example  of  the  analysis  needed  to  justify  (30)  is  given  in  a  later  para- 
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graph  for^the  two-pulse  staggered  case.  Assuming  that  this  approximation 


is  valid  then  the  average  power  spectrum  of  the  clutter  process  can  be 
written  as: 
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(31) 

where 


(32) 


denotes  the  average  clutter  power  per  range  ring.  It  is  shown  in  the 
Appendix  that 


M-l 
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k,0 


(33) 


hence,  the  clutter  spectrum  reduces  to 


Sc(f) 


^•77-  i 

s  p  i°-0 


vf  -  b 


(34) 


Receiver  Noise  Model 

It  follows  directly  from  I-(36)  that  staggering  the  transmitter  PRF 
has  no  effect  on  the  receiver  noise  process.  Therefore,  it  remains  a  zero- 
mean  white  noise  process  with  spectral  dersity  2NQ. 


Two-Pulse  Staggering 


In  order  to  gain  some  physical  understanding  of  the  mathematical 
expressions  for  the  target  and  clutter  spectra  the  special  case  of  a  two-pulse 
stagger  will  be  studied.  This  is  illustrated  in  Figure  2.  Using  M  =  2, 
eo  =  0»  ei  =  e  08)  and  (19)  the  transform  of  the  target  signal  is 


S(f;a)  =  y  e 


*j2tff  T 


hi 

P  i=-o 


Cg(x»v)Fg(f  -v  -  y~)  +  C-j  (t,v)Fg(f-vf^| —  y-) 


(35) 


where 


Cq(x,v) 

C^xtv) 


~  j  2ttv£ 

O  +  e  )/ 2 


|I(T) 


-J2ttv£ 
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(36a) 

(36b) 


Hence  the  spectrum  of  the  target  return  is 


(cos^irve)  jF_(f-\r  t~)|^  ♦  (sinSrvc) 
1  y  P 


Fg(f-v 
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V 


(37) 


Typical  plots  of  the  target  spectrum  are  illustrated  in  Figure  3.  There  are 
two  significant  observations  to  be  made:  (1)  whereas  in  the  uniformly  sampled 
case  all  of  target  energy  is  located  at  PRF  multiples  of  the  true  Doppler* 
staggering  causes  the  energy  to  be  split  into  two  pieces  separated  by  one-half 
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the  PRF,  and  the  pair  folds  over  at  the  PRF,  and  (2)  whereas  in  the  uniformly 
sampled  case  targets  moving  at  dopplers  greater  than  a  PRF  led  to  spectra  that 
were  indistinguishable,  now  the  fundamental  ambiguity  occurs  with  a  period 
2 ft.  This  shows  that  staggered  PRF's  provide  a  basis  for  unambiguous  velocity 
estimation. 


From  (31),  the  exact  fora  of  the  clutter  spectrum  reduces  to 


sc(f) 


a2 


V 


I 

COS^tr{fr-  y-)c 

|F9!f 

'V 

i  =  .ao 

P 

Since  the  frequency  extent  of  F  (f)  is  very  narrow  relative  to  1/T  ,  it  can 

V  r 


reasonably  be  assumed  that 

(cos2-fe)  fg(f)  2  .  |fg(f)  2 

(39a) 

(sln2»fc)|Fg(f)  2  •  (sfsl^tFjj2 

(39b) 

Using  these  approximations  the  clutter  spectra  can  be  sketched  as  shown  in 
Figure  4  from  which  it  is  observed  that  as  for  the  target  spectrum  the  clutter 
power  also  splits  into  two  pieces,  one  piece  being  located  at  DC,  the  other 
at  -  1/2  Tp,  with  the  aggregate  folding  over  at  the  PRF.  The  simple  sketch 


has  been  drawn  to  indicate  that  the  clutter  power  at  1/2  Tp  is  significantly 
smaller  than  that  at  OC.  A  quantitative  measure  of  the  relative  power  in  each 
of  these  terms  can  be  found  by  integrating  (39a)  and  (39b).  This  has  been 
done  for  the  sin  x/x  antenna  pattern  and  ARSR  system  parameters  and  it  was 
found  that  the  clutter  power  at  1/2  Tp  56  dB  down  from  that  at  DC.  Since 
10  bit  A/D  converters  correspond  to  a  subclutter  visibility  no  greater  than 
4G  dB,  the  effect  of  the  clutter  power  at  1/2  T  is  negligible,  hence  justi- 

r 

fying  the  assumptions  leading  to  the  clutter  spectrum  in  (34). 

Therefore,  '.ne  implications  of  the  staggered  PRF  are  now  clear:  Whereas 
the  energy  of  a  moving  target  return  splits  into  two  pieces,  the  clutter  power 
continues  to  fold  over  at  multiples  of  1/Tp.  Hence,  if  the  target  Doppler  is 
also  a  miHiple  of  1/Tp,  namely  a  former  blind  speed,  then  although  one  por¬ 
tion  of  the  target  energy  is  masked  by  the  DC  clutter,  the  other  portion  is 
located  in  a  relatively  clutter-free  area  at  1/2  Tp.  This  is  the  reason 
staggered  PRF  enhances  target  detection.  However,  if  in  addition,  filters 
that  are  matched  to  the  target  spectrum  are  constructed,  then  it  appears 
that  Doppler  estimation  over  a  frequency  interval  larger  than  one  PRF  is 
possible.  Although  the  topic  is  discussed  in  more  detail  in  Section  IV  we 
briefly  discuss  the  implementation  of  the  filter  matched  to  the  two-pulse 
staggered  signal  spectrum. 

Matched  Filter  Realization 

From  the  preceding  discussion  it  is  shown  that  the  target  spectrum  is 
a  unique  function  of  the  true  target  Doppler  over  an  Interval  that  can  be 
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many  times  larger  than  the  PRF.  If  a  matched  filter  bank  could  be  construe 
ted  then  not  only  would  the  detection  performance  be  optimized  but  unambigu 
ous  estimation  of  target  Doppler  would  be  possible. 


Let  us  suppose  that  we  require  the  resolution  of  velocity  to  within 
the  interval  Av  =  1/NT  .  Then  each  PRF  interval  can  be  quantized  into  N 

r 

subintervals  and  we  can  then  express  the  true  target  Doppler  as 


_  ->0  +  no 

v  “tt  nt  • 

p  p 


From  (35)  the  signal  spectrum  for  the  two-pulse  stagger  is 


=  y  X  co(T-v)Fg(f-  sr  -  -t r) 

•  r  r 

1  =  ^oo 


+C1  (t>v)F  (f-  -  ~J~~  +  2T~ ) 

y  p  p  i» 


The  infinite  sum  shows  the  periodic  foldover  of  the  target  spectrum  at  the 
PRF.  From  a  measurements  point  of  view,  nature  allows  us  to  observe  this 
function  only  in  the  interval  [0,  1/T  ].  Then  we  see  only  the  function 


S(f;a)  =  y  Cq(t,v)F  (f-  +  C^(t,v)F  (f- 

P  P  P 


We  can  readily  construct  a  bank  of  filters  that  extend  over  the  [0,  1/Tp] 
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range  where  each  filter  is  tuned  to  the  function  F  (f-  n=0,1,...,N-l. 

By  themselves,  these  are  not  matched  to  the  specified  signal.  To  accomplish 
this,  we  combine  weighted  pairs  of  filters  that  are  separated  by  l/2Tp  Hz. 
For  the  filters  tuned  to  n/NT  and  (n-l/2}/NT  we  apply  the  weights 

r  r 

C*(t,  y~  +  ,  C*(r,j-  +  for  i=0,+l ,+2,...,  +M.  For  each  value  of 

P  P  P  P 

i,  this  gives  rise  to  another  filter  with  transfer  function 


iao,+,...,  +M,  n=0,l,...,N-l 


When  i=io,  nanQ,  this  filter  is  matched  to  the  two-pulse  staggered  signal. 


F  (f-  -L-) 

,Vf  NTp; 


From  a  practical  point  of  view,  the  sub-bank  of  filters 
N-l 

can  be  formed  by  taking  an  N-point  Discrete  Fourier  Trans - 


n=0 


form  (DFT)  of  the  received  signal.  The  super-bank  of  filters  is  then  ob- 

1L  ^  J 

tained  by  multiplying  the  ntn  DFT  coefficient  by  Cq(t;  and  the 


(n-  L~T  coefficient  by  C*(t;  +  jjy-)  for  1=0,  +1,...,  +M.  Therefore, 

P  P 

an  N-point  TF”  r;ives  rise  t,  a  bank  of  2MN  matched  filters  that  extend  over 


M  M 

the  frequency  Interval  |-  y~*t 

P  PJ 

coefficients  in  the  right  way. 


simply  by  combining  the  outputs  of  the  DFT 
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In  Section  IV  we  return  to  this  discussion  in  more  detail  when  we 
consider  the  MTI  ambiguity  function.  In  Section  III, a  quantitative  measure 
of  the  improvement  in  detection  performance  will  be  evaluated  using  the 
Signal-to-Interference  Ratio  (SIR)  that  was  derived  in  Part  I. 
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III.  SIR  PERFORMANCE  ANALYSIS  FOR  STAGGERED  PRF 


The  SIR  for  an  arbitrary  linear,  sarapled-data  filter  when  sampled  at 
time  t  was  givenby  I-(72),  viz. 

1/2T 


,  ,P  j2irfT  j2irfT  j2Trfx 

/  H(e  p)Zs(e  p;2o)(e  ) 

'-1/2T 

p(t)  =  Tp  - 


df 


2 


/  V^V 


-1/2T 


j2irfT0 

Sc(e  p)  ♦  2NQ 


df 


(40) 


Even  though  the  transmitter  PRF  is  staggered,  the  sampled-data  processor 
operates  on  the  samples  of  the  signal  and  noise  and  it  natters  not  when  those 
samples  were  taken.  Therefore,  (40)  applies  to  the  present  problem,  although 
it  is  noted  that  the  signal  spectrum  will  be  different,  due  to  the  non-uniform 
sampling.  As  before,  it  is  noted  that  only  those  frequency  terms  in  the  in¬ 
terval  (-l/2Tp,  l/2Tp)  are  of  interest.  This  is  consistent  with  (19)  since 

the  frequency  dependence  shows  up  only  in  terms  like  F_(f-V  ♦  -  i-) 

9  'p 

which  is  folded  over  every  1/Tp  hz. 


Using  the  Schwarz  inequality  it  is  easy  to  show  that  (40)  is  maximized 
by  choosing 

z  *.  j2*fTr 

’>  *  — ' “TCT - —  (41) 


j2*fT  l  (e 
HCe  p'  -  s 


WT 

Sc(e  p)  +  2Nc 
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which  is  the  clutter  filter,  matched  filter  cascade  combination.  When  this 
Is  done,  the  resulting  maximum  value  of  the  SIR  is 


popt  “  Tp 


P  Z$(e 


1/2T 

1  SCe 

•l/2Tp  Ve 


j2irfT 


TSSfT 


df 


P)  +  2N 


(42) 


The  aliased  clutter  spectrum  is  given  by  (34),  but  since  the  integration  ex¬ 
tends  over  the  (-  1/2Tp,  1/2Tp)  frequency  interval,  only  the  term  about  DC 
need  be  taken.  Taking  the  squared  magnitude  of  (19)  and  using  the  approxlma 
tion  in  (27),  the  target  spectral  density  reduces  to 


(43) 


Using  these  results  and  the  fact  that 


ek(x,v) 


ak(?)bk(v) 


(44) 


which  follows  from  (18),  then  the  SIR  in  (42)  becomes 


H-l 

‘>opt(vo)  ‘  Kr  1  ImvI*  /  "  T2-: 

k*g  0  ’ 


,/2Ip  ? 


Fg(f  "vo  *  *  T* 

P  P 


df 


-1/2tptt  v  ivf)r 


(45) 
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Rather  than  attempt  a  rigorous  evaluation  of  (45) r  it  is  easier  to  draw  upon 
the  physical  understanding  of  the  target  and  clutter  spectra  to  simplify  the 
SIR  expression.  It  was  shown  in  the  last  section  that  the  M-pulse  staggered 
PRF  causes  the  target  energy  to  split  into  M  components  that  are  folded  over 
into  the  (-l/2Tp,  l/2Tp)  interval,  while  the  clutter  was  distributed  about 
DC.  Since  the  frequency  extent  of  Fg(f)  is  narrow  relative  to  the  window 
l/MTp,  there  are  values  of  for  which  there  is  no  interaction  between  the 
clutter  and  target  spectra.  In  this  case,  for  each  vQ  there  is  a  value  of  i 
that  puts  Fg(f  -  vQ  —  |-)  within  the  (-l/2Tp,  l/2Tp)  interval  and 


df 
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In  the  Appendix  it  is  shown  that 


n-  i 

|Ibk(V 


whence  it  follows  that 


popt  (V  - 


Y  ! 2  E 

70|  5 


This  is,  of  course,  just  the  coherent  integration  gain  provided  by  matched 
filtering  the  target  out  of  the  white  noise  background. 

The  SIR  degrades  from  this  optimum  value  when  any  one  of  the  H 
components  of  the  target  spectrum  interacts  with  the  clutter  spectra.  The 
worst  case  occurs  when,  for  some  k  and  i,  kQ  and  iQ  say, 

k  1 

-y0  *  ht  -  r 3  0  (so) 

p  p 


In  this  case,  since  the  dutter-to-white  noise  ratio  is  very  large. 


7Te  K(f  '  vo  '^’l2  Ts 

I  K - E —  a 

-lmp  f  Fg(0  2  ♦  2N0Tp  0 


For  the  remaining  H-l  components  of  the  target  spectrum  that  are  located 
within  the  ( -l/2Tp,  l/2Tp)  interval,  there  Is  little  interaction  with  the 
clutter  spectra.  Hence,  for  those  values  of  k  t  k  (46)  holds  and  the  SIR 
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can  be  written  as 


popt(V 


c  M-l 

2nV  £ 

0  P  k=0 

m. 


bk(v0) 


2  +  ?k„('’°) 

0  o 
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2e. 


0  P 


2NoVs 


o  g 


2r 


*rr 

0  p 


1- 


bkft {vo} 
0 


(52) 


where  the  last  approximation  follows  from  the  fact  that  the  clutter  to 
receiver  noise  ratio  is  »  1.  This  expression  for  the  SIR  holds  for  values 
of  VQ  given  by 


K 

o.o 

sr  r 


(53) 


where  first  a  value  of  mQ  is  chosen  and  then  for  each  m0,  kQ  =  0,1 ,2»* * ♦ ,M-1 . 
Then  the  optimum  SIR  performance  curve  can  be  sketched  by  using  the  formula 
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(55a) 


For  the  case  of  a  two-pulse 

IW2 

IW2 

so  that 


r" 

2 

^0 

E9 

stagger,  M  =  2,  kQ  =  0  or  1  and 

=  C0S2TnjQ£ 

2 

=  sin  ttuq£ 


m 

o 


2 

sin  irue 

0 

if 

< 

o 

u 

m 

cosVoQe 

if 

0 

vo=Tp 

1  otherwise 


(55b) 


(56) 


Whereas  when  no  pulse  staggering  is  used  (e  »  0),  the  SIR  is  essentially 
zero  at  multiples  of  the  PRF,  staggered  pulse  transmissions  lead  to  mean¬ 
ingful  detection  performance,  especially  at  higher  Doppler  velocities.  The 
price  paid  for  this  enhanced  performance  at  the  blind  speeds  is  a  degradation 
in  the  SIR  performance  at  intermediate  Doppler  frequencies.  These  results 
are  summarized  in  the  SIR  performance  curve  plotted  in  figure  5.  It  is 
worth  noting  that  similar  results  can  be  obtained  for  the  pulse  canceller 
clutter  filters  by  working  directly  from  (40)  using  the  appropriate  filter 
transfer  functions.  The  SIR  performance  of  the  ASR-7  that  uses  a  6-pulse 
stagger  algorithm  is  shown  in  Figure  6. 
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performance  usirvg  an  operational  tix-puise  stagger. 


IV.  STAGGERED  PRF  AMBIGUITY  FUNCTION 

In  Section  II,  the  target  spectrum  resulting  from  a  staggered  PRF 
transmission  sequence  was  derived  and,  for  the  two-pulse  case,  illustrated  in 
Figure  3.  It  was  noted  that  the  spectra  for  targets  separated  by  Doppler  shifts 
greater  than  one  PRF  were  not  identical  as  was  the  case  when  uniform  sampling  was 
used.  This  indicates  that  it  may  well  be  possible  to  estimate  target  Doppler 
unambiguously.  This  question  is  most  easily  examined  by  evaluating  the 
ambiguity  function  of  the  staggered  PRF  pulse  train.  The  calculation  is  not 
conceptually  difficult  but  it  can  become  tedious.  In  order  to  develop  some 
intuition,the  clutter-free  ambiguity  function  will  be  computed  first  and  then 
generalized  to  the  situation  in  which  the  clutter  filter  is  present.  In  the 
former  case  the  ambiguity  function  is  the  delay-Doppler  distribution  of  the 
output  of  the  matched  filter.  It  is  denoted  by  where 

1/2  Tp 

*  ./  2S  Cfi&)  df 

-1/2  T  (57) 

p 

Rather  than  attempt  to  evaluate  (57)  by  direct  substitution  it  is  easier  and 
core  instructive  to  draw  heavily  upon  the  physical  interpretation  of  the 
correlation  operation  implied  by  this  equation.  The  necessary  intuition  can 
be  developed  by  studying  the  transmitted  signal  for  the  three-pulse  staggered 
case.  From  (18)  and  (19)  the  transform  of  the  target  signal  is 
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(58) 


The  magnitude  characteristic  of  this  function  is  illustrated  in  Figure  7a. 

It  will  be  assumed  that  t,  tq  and  vQ  are  fixed  so  that  the  correlation  opera¬ 
tion  in  (57)  can  be  studied  as  a  function  of  v.  Kaking  use  of  the  1/Tp 
periodicity  in  the  target  spectrum,  the  integral  in  (57)  can  be  evaluated 
using 


h  {f;  T’v)  Zs  (f;  V  vo5df 


v  -1/T 
o  p 


(59) 


The  first  situation  of  interest  occurs  when  v  =  v  in  which  case  the  Ooppler 

o 

coefficients  line  up  exactly.  Equation  (59)  becomes 
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(60) 


Assuming  the  t  takes  on  only  integral  values  of  Tp, 


then  from  (18a) 


ak(x)  =  e 


(61) 


It  then  follows  that 
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(62) 


ittYn~r^v  ■ aaaaaiaiitiBi^mii^^  . - . . . 


where 


1/2T, 
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j  2ufT 


df 


(63) 


-1/2T. 


In  Part  I  it  was  shown  that  this  function  was  precisely  the  autocorrelation 
function  of  the  two-way  antenna  pattern.  Then  the  ambiguity  function  when 
the  signals  are  matched  in  Doppler  is 


^•vvV 
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where  from  (18b) 
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It  is  shown  in  the  Appendix  that 
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hence 
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a  result  which  Is  Intuitively  satisfying. 


In  order  to  evaluate  the  ambiguity  function  at  other  values  of  v»  It 
Is  useful  to  think  of  gradually  increasing  v  from  Its  value  at  vQ.  For 
example,  when  vQ  <  v  <  1/3  Tp,  the  absolute  value  of  Zs(f;T,v)  Is  shown  In 
Figure  7b.  When  the  correlation  operation  Is  performed  to  evaluate  (59) 
for  these  values,  there  will  be  no  spectral  overlap  and  the  ambiguity  function 
will  essentially  be  zero.  No  significant  contribution  will  be  made  to  the 
ambiguity  function  until  v  *  vQ  +  l/3Tp.  In  this  case,  different  frequency 
coefficients  line  up  and  (59)  becomes 
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From  (61)  it  follows  that 
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and  therefore 
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The  ambiguity  function  Is  then 
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The  next  step  is  to  set  v  °  v  +  xy-  and  repeat  the  above  operations.  In  this 

o  j  ip 

case,  the  coefficients  are  displaced  by  two  and  the  ambiguity  function  becomes 
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This  process  continues  ad  infinitum  and  it  is  possible  to  deduce  a  rule  for 
generating  the  ambiguity  function.  In  the  M-pulse  stagger  case,  it  becomes 


n- 1 

^T*To*vo  +  S5r*vo)|  53  ^4Rg(T'To}|2  WV  “)bk(V 


where  for  positive  values  of  n  =  0,1,2,’*’,  m  takes  on  the  values 
m  «=  0,1 ,2,’ ” ,M-1 ,  and  for  negative  values  of  n  e  m  «  M-l,  M-2,  0. 

In  (73)  use  has  been  made  of  the  fact  that 


Wv)  a  b(k«)  M 

modulo  m 


which  follows  directly  from  (65).  It  is  shown  in  the  Appendix  that  (73)  can  be 
reduced  to 
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,  m+nM  \ 
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p  s> 


k=o 


(75) 


which  is  a  function  only  of  the  differences,  t-tq  and  v~ ;Q,  and  the  stagger 
parameters  eQ,  e-| ,  ***,  eM_i*  It  can  be  deduced  immediately  that  the  ambiguity 
function  is  unchanged  if  eQ  =  o,  hence  for  an  M-pulse  stagger  there  are  M-l 
parameters  that  can  be  chosen  to  shape  the  ambiguity  surface. 


For  the  special  case  of  a  two-pulse  stacger  (75)  reduces  to 


C(t-T0,V-V0) 
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sin  it 
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V*p 


(76) 


and  this  is  sketched  in  Figure  8a  and  compared  with  the  ambiguity  function 
for  the  uniformly  sampled  case  in  which  t  a  o,  in  Figure  8b.  It  is  clear 
therefore  that  Doppler  i-esolutlon  i-.  theoretically  possible.  Whether  or 
not  the  stagger  parameters  can  be  chosen  to  force  the  subsidiary  side-lobes 
below  a  practically  useful  level  is,  however,  a  separate  question.  It  is  of 
interest  to  examine  the  ambiguity  function  of  higher  order  stagger  sequences 
that  are  currently  used  in  practice.  The  results  for  the  ASR-7  radar,  that 
uses  a  6-pulse  stagger  are  shown  in  Figure  9. 
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(a)  Two-pulse  staggered  PRF. 
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V  -  vn 


(b)Unifomi  PRF. 


Fig.  8.  MTI  ambiguity  function. 


The  Matched  Filter-Clutter  Filter  Ambiguity  Function 


In  the  preceding  section,  the  ambiguity  function  for  the  clutter-free 
case  was  derived.  This  is  a  useful  characterization  when  the  signal  is 
designed  to  function  in  only  a  white-noise  environment  as  it  is  then  clear 
that  all  of  the  side-lobes  should  be  made  uniformly  low.  The  more  typical 
situation  for  MTI  requires  a  characterization  that  includes  the  clutter  in 
the  analysis.  If  the  ambiguity  function  is  viewed  as  the  delay-Doppler 
energy  distribution  of  the  signal  out  of  the  optimum  processor,  then  it  is 
clear  that  the  effect  of  the  clutter  is  to  add  notch  filters  at  multiples 
of  the  PRF's.  Then  the  more  general  ambiguity  function  is  given  by 


1/2TP 

/  Hc(f)zs  (f;o)Z^figo)df 


(77) 


As  in  the  clutter-free  problem  the  general  result  will  be  obtained  by  extending 
the  arguments  made  for  the  three-pulse  staggered  case.  This  is  most  easily 
done  by  writing  a  general  expression  for  (60)  and  (68)  from  which  the 
ambiguity  function  is  deduced.  This  expression  is 
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(78) 


For  the  purpose  of  this  discussion  it  is  reasonable  to  assume  that  the 
clutter  filter  transfer  function  changes  slowly  over  the  width  of  the  signal 
spectrum,  hence  allowing  the  following  approximation  for  the  last  term  in 
(78) 


where  the  last  equation  follows  from  a  generalization  of  (60)  and  (68).  Then 
the  ambiguity  function  is 
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(80) 


To  evaluate  (80)  it  is  assumed  that  the  clutter  filter  is  well  modelled  by  a 
notch  at  DC  as  well  as  at  all  multiples  of  the  PRF.  The  approximation  was 
developed  in  conjunction  with  the  evaluation  of  the  SIR  for  staggered  PRF's. 
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Suppose  now  that 


v0  "  ffy~ J*  j~  k  3  0,1  ."-.M-l  (81) 

P  P 

holds  for  every  n,  then 

Mvc  “  HlM  “  1  (82) 

and  (80)  reduces  to  the  clutter-free  ambiguity  function.  If  for  some  value 
of  k,  k*  say, 


*  sir 3  f~  {83) 

p  p 

for  some  value  of  n,  then 


"c^o  -  '  0 


(84) 


and  (80)  reduces  to 


This  function  is  much  more  difficult  to  plot  as  it  depends  on  the  true  target 
Doppler  rather  than  just  the  difference  between  the  true  and  tested  values. 
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In  fact,  for  an  M-pulse  stagger,  (W+l )  cuts  of  the  ambiguity  function  are 
needed  to  describe  it  completely. 


MTI  Signal  Design 


In  the  clutter-free  case,  it  is  clear  that  the  stagger  parameters 
should  be  chosen  tc  produce  an  ambiguity  surface  ;dth  uniformly  low  Doppler 
side-1  rtbes.  From  (73)  this  reduces  to  the  problem  of  picking  the  stagger 
parameters  eQl  ^ so  that 


^M-l  *ck)f 

|r2  e  p  | 


(86) 


where  fcr  ZGdB  sidelobes  6  would  be  *1,  etc.  This  signal  design  problem 
has  fitoch  in  coeason  with  design  of  antenna  patterns  using  an  array  with  non- 
uniforaly  spaced  elements.  This  is  a  difficult  problem  to  solve  and  it  is 
expected  that  when  the  clutter  filter  is  added,  it  would  be  even  more  difficult 
to  simultaneously  design  the  (K+l)  folds  of  the  cluttered  ambiguity  function. 

A  simpler  design  strategy  can  be  obtained  from  the  SIR  analysis  in  Section  III 
where  it  was  shown  that  the  degradation  in  the  performance  was  given  by  (50). 
From  this  expression  it  is  clear  that  ihe  stagger  parameters  could  be  chosen 
to  minimize  the  depths  of  the  notches  by  minimizing  the  functions 
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which  follows  from  the  definition  of  bm(v)  in  (65).  This  expression  can 
easily  be  manipulated  to  take  exactly  the  same  form  as  that  in  (86).  This  is 
interesting  as  it  shows  that  the  simple  criterion  of  uniformly  low  sidelobes 
is  a  good  signal  design  strategy  in  the  cluttered  as  well  as  the  white  noise 
environments. 


Unfortunately,  time  did  not  permit  the  thorough  examination  of  these 
signal  design  problems.  Therefore,  as  of  this  writing,  their  solution 
remains  an  unanswered  question  and  it  will  be  necessary  to  be  content  to  use 
the  MTI  ambiguity  function  as  a  tool  for  signal  analysis,  and  only  indirectly 
for  signal  synthesis. 
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APPENDIX 


In  order  to  evaluate  the  MTI  ambiguity  function  it  was  claimed  in 
(75)  that 
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where  from  (65) 
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It  is  the  purpose  of  this  section  to  prove  (A-l)  and  hence  deduce  that 
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kao 
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(A-3) 


as  stated  in  (48  and  (66).  To  begin,  vectors  and  matrices  are  defined  as  follows: 


/  ~j2irve0  -J  2*061  ”j2rrveM  ,, 

&(v)  =  (e  e  e  M-^ 


(A-4) 


m,  x 
X  M  13 


ifl+M-1 


(A-5) 


48 


-j2ir{k+m)j=j- 

(«")kp  =  ^  (A-6) 

where  m,  k  and  p  take  on  values  0,1 ,2, *  *  * ,M-1 .  First,  from  (A-2)  it  is  noted 
that 
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where  use  has  been  made  of  the  fact  that 
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Hence 
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If  <u_, 


where 


Then 


v>  denotes  the  inner  product  in  C^,  then 


,  /  ,  m+nM  s  _  ,  m,  ,  m+nM  %  o,  , 

*(vo  +  "mT  *V  "  <  *  (vo  +irk  W 


=  <  M^Vq  +  M°B(v0) 

=  <  M°  ^(vq  +  ,£(v0)  > 

P 


denotes  conjugate  transpose.  Now  let 
Qni  =  M°*Mm 


(A-10) 


(A-1 1 ) 


PA 


3 


p=o 


50 


(A-12) 
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For  the  case  when  k^j,  the  second  term  in  (A-12)  can  be  evaluated  by  setting 
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(A-l 4 ) 


and  noting  that 
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Hence  =0  when  kjty  and  therefore  the  matrix  Qm  is  diagonal  for  all  m. 


Using  this  result  in  (A-l 0)  yields 
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as  required. 
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